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Abstract 

Using the optical reference geometry approach, we have derived in the fol- 
lowing, a general expression for the ellipticity of a slowly rotating fluid con- 
figuration using Newtonian force balance equation in the conformally pro- 
jected absolute 3-space, in the realm of general relativity. Further with the 
help of Hartle-Thorne (H-T) metric for a slowly rotating compact object, we 
have evaluated the centrifugal force acting on a fluid element and also evalu- 
ated the ellipticity and found that the centrifugal reversal occurs at around 
R/R s ~ 1.45, and the ellipticity maximum at around R/R s ~ 2.75. The 
result has been compared with that of Chandrasekhar and Miller which was 
obtained in the full 4-spacetime formalism. 

1 Introduction 

One of the enigmatic problems in the context of pulsars is still the under- 
standing of the internal structure of rotating compact objects. As normally 
one considers the object to be slowly rotating, most of the model calcula- 
tions have considered the well known approximate solution of Hartle and 
Thorne Jl| , |J as the basic solution for the fluid configuration in the context 
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of general relativity. It is generally believed that for any such rotating fluid 
configuration if one considers the 'force balance', in the purely Newtonian 
physics one encounters no strange behaviour irrespective of the size of the 
compact object as the two traditional rivals the gravitation vis a vis the cen- 
trifugal force acting on a fluid element would always be opposing each other. 
However, in the context of general relativity as has been shown by Abramow- 
icz and Prasanna || for a sufficiently small size object, the centrifugal force 
acting on a test particle of mass m in circular orbit outside a static mass M 
is given by 

where v is the speed of the particle as seen in the conformally projected 3- 
space of the optical reference geometry (ACL) [|J. As seen from the above 
expression the centrifugal force would not oppose gravity if the particle is 
situated at a distance r < 3M. As there could exist ultra compact bodies 
[§] of size 2M < R < 3M, it would become relevant to consider the effect of 
such a centrifugal force reversal on a fluid element of a possible ultra compact 
rotating configuration. 

Another important manifestation of the same result viz, introducing New- 
tonian forces in general relativity is the explanation of ellipticity maximum 
for a rotating configuration given by Abramowicz and Miller ||, an effect 
discovered by Chandrasekhar and Miller Though the explanation given 
by Abramowicz and Miller is qualitatively viable, quantitatively there ap- 
pears a difference in the location of the ellipticity maximum, which perhaps 
is due to the fact that they considered only the Schwarzschild background 
geometry, which does not take into account the influence of rotation of the 
central body inherently. 

In the following, we reexamine the scenario by studying the possible cen- 
trifugal reversal and the ensuing ellipticity maximum for a slowly rotating 
fluid configuration by adopting the Hartle - Thorne solution which indeed is 
better suited to study slowly rotating fluid configuration. 
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We start with a general axisymmetric, stationary fluid configuration and in- 
troduce a formalism to treat the four forces on a fluid element in the 3+1 
conformal splitting and then adopting Hartle's solution as a specific example 
consider the centrifugal force. Using the Newtonian principle of force balance 
equation for a rotating spheroid we then derive a general expression for the 
ellipticity and again study its behaviour for the Hartle solution. We find 
that the result matches closer to that of Chandrasekhar and Miller result 
thus validating the more general expression derived. 



2 Formalism 

The equation of motion for a perfect fluid distribution on a general curved 
space time 

ds 2 = gijdx % dx\ (2) 

are given by 

(/> • /'Hr: ; rn h ij P>j (3) 

where p is the matter density, p the pressure, U l the four velocity and 
the projection tensor (g 1 ^ + U l U^). This may indeed be expressed as the four 
force acting on a fluid element 

fi := (,> ■ i>)ir, :j r j ) ■ h>,,, (4) 

= {p + p)[U^d j U i -^U m U^d i g jm ]+h? i p, j , (5) 

which when p = and p = mo, reduces to the well known four force expres- 
sion acting on a particle || 

mofi = P^P, - ^P m d i9jm . (6) 

Using the ACL formalism of 3+1 conformal slicing of the space time 

ds 2 = dl 2 - g 00 (dt + 2cu a dx a ) 2 , (7) 
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with dl 2 representing the positive definite metric of the absolute 3-space 
g^dx^dx", equation (5) may be rewritten as 



fa 



fo = ^- 1 (p + p)U^U + h^ 
+2u a f + (h£ - 2u a h$)p tli 



(8) 
(9) 



with 



M 2 






= W 




= 9<jjt 






goo 


= -$ 


g^o 


= -2$^ 


g^v 





(10) 



and V M denotes the covariant derivative in the absolute 3-space. 

As we are interested in slowly rotating fluid configurations, we shall consider 
the most general stationary axisymmetric space time metric as given by M 



(11) 



where u, ip, fix, and u are functions of r and 9 and dcf) = d<p — udt. For this 
metric when we have only = ■£ ^ 0, the normalisation condition yields 



.-2v 



-1/2 



(12) 



where Q 



iii 



is the angular velocity of the fluid as measured by the station- 



ary observer, and uo is the angular velocity acquired by an observer falling 
freely from infinity As this frame (t, <p, r, 8) is static, the equations (8) and 
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(9) simplify considerably and then one can calculate the 3-component of the 
force as given by 



/« = $ (P + P) 



M 2 



(13) 



which in fact, when zero, gives the equation of hydrodynamic equilibrium for 
a rotating fluid configuration. Now given a metric (approximate or post- 
Newtonian solution) one can calculate the 'centrifugal acceleration' term 



M 

£7 M V 'fj,U a and 'gravitational acceleration' term -^-d a ^, which for the met- 
ric (11) yields 



cf 



e 2 V, 



(14) 



where prime denotes differentiation with respect to r. 
Ellipticity 

It is well known that a rotating fluid configuration deviates from spherical 
symmetry and depending upon the degree of rotation the equatorial diameter 
expands whereas the polar diameter contracts thereby producing a change 
in shape. The various equilibrium configurations of rotating fluids have been 
well discussed in the literature, and the sequence goes through Maclaurin 
spheroids to Jacobi ellipsoids [[J. For slowly rotating configuration one can 
consider the Maclaurin spheroid with the ellipticity defined through the usual 
definition of 



- (1 - e 2 ) 1 / 2 

l_ e 2)l/6 : 



(15) 



e being the eccentricity defined as e 
and equatorial radii respectively. 



'1 —tf/a 2 ) 1 ' 2 , where b and a are polar 
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Maclaurin had shown that the acceleration due to gravity at the equator and 
pole has the values 

9eguator = 2itGpa ^p^ [sin~ l e - e(l - e 2 )^] 

(16) 

= ATrGpa^P^-ie-il-^sin^e], 



9poi 



e 



wherein he had considered the possible effects that could arise due to the 
internal stresses in the body. However, as we are looking for a solution in 
general relativity, wherein the gravitational field inside the body is described 
through a metric which is a solution of Einstein's equations for a perfect 
fluid distribution, the gravitational potentials gij would be incorporating the 
effects of all characteristics of the fluid distribution . With this proviso, in 
the new language of optical reference geometry it is sufficient to consider 
the modified expression for the gravitational and centrifugal accelerations 
as given by (14) and use the Newtonian force balance equation to relate the 
eccentricity with the acceleration. Thus generalising the Newtonian equation 

9 equator ~ ^ = 9pole ( 1 - (? ) * • (17) 

to 

F ge -F cf = F gp (l-e 2 ) 1/2 (18) 
and using the force expression as given by 

(9 = n/2):F ge = e^^/ 2 V (r, vr/2), (19) 

(6 = 0):F gp = e 2 ^°^(r,0), (20) 

and F c f as in (14), the eccentricity of the configuration would be given by 




^cf Fge 

F 



(21) 



and the ellipticity in the limit of slow rotation e << 1, 
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Hartle's Solution 



For a slowly rotating perfect fluid configuration Hartle has obtained an ap- 
proximate solution of Einstein's equations given as 



ds 2 = [-e»°(l + 2(h + h 2 P 2 ))]dt 2 +(l 



2M 



' (mp + m 2 P 2 ) 
1 + 2 (r-2M) 



dr 2 + r 2 [l + 2(v 2 -h 2 )P 2 ] ( 23 ) 
dfl 2 + sin 2 6{d<f) - udt) 2 ] + O (Q 3 ) , 

which represents the rotation as perturbation over the non-rotating metric 

ds 2 = -e"°(r)dt 2 + [1 - 2M{r)]~ 1 dr 2 + r 2 [dQ 2 + sin 2 Bd$ 2 ). (24) 

Expressing the forces as obtained above in terms of the H-T potentials one 

gets 



F cf = r 2 u 2 (l/r-u' /2), 



Fge = -e" [v' (l + 2h -h 2 ) + 2ti -ti 2 }, 



F 9P = -e U0 [u' (l + 2h + 2h 2 )+2h' + 2h' 2 }, 



yielding for the ellipticity 



r 2 u 2 



e = 3(h 2 + h' 2 /u' ) + — (2/ru' - 1) . 



(25) 
(26) 

(27) 
(28) 



3 Results and discussion 

In the present work we have derived the general expressions for the centrifugal 
acceleration (F c f) (using ACL formalism) and ellipticity (e) (replacing New- 
tonian acceleration by relativistic counterparts in the force balance equation 
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(17)). Using H-T solution we obtain the values of the centrifugal force and 
the ellipticity for a slowly rotating configuration. For the sake of comparison 
we have also calculated the ellipticity as given by Chandrasekhar and Miller 



( \ 3 



(29) 



where 



6 = 2p* 2 /(dv /dr) , p* 2 = -h 2 - l -r 2 e~^u 2 . (30) 



Writing the expressions in dimensionless units: 



cf 



(G 2 J 2 /c 4 R 5 s ) ' (G 2 J 2 /c 6 Ri)' 



(31) 



where J is the angular momentum and R s (= 2M) is the Schwarzschild radius, 
we have evaluated the quantities for a series of homogeneous configurations 
with decreasing radii keeping M and J conserved, e denotes our calculations 
(equation (28)) whereas Ih-t denotes the values for corresponding configu- 
rations as obtained using equation (29) in table 1. 

Comparing our present result with that of Abramowicz and Miller, who had 
obtained the maximum at R/R s = 3, using pure Schwarzschild geometry, we 
see that incorporating the effects of rotation in the geometry (even approx- 
imately) improves the result as the maximum R/R s ~ 2.75 shifts closer to 
that obtained by Chandrasekhar and Miller R/R s ~ 2.3 (fig 2). 

Regarding the centrifugal force the general expression obtained above does 
show the reversal at R/R s ~ 1-45 and a maximum at R/R s ~ 2.1 (fig 1). It is 
interesting to note that even after including the effects of fluid distribution in 
the space time geometry, the centrifugal force reversal seems to occur closer 
to the value as was known in the Schwarzschild geometry. However, as the 
ellipticity maximum indicates a possible change in shape of the configuration, 
it is to be noted that our expression shows that for a collapsing configuration, 
the change occurs earlier (R/R s ^ 2.75) than what had been obtained by 
Chandrasekhar and Miller (R/R s w 2.3). As the shape of the body does 
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depend upon the ellipticity as its value starts decreasing after reaching a 
maximum, the body would in principle tend towards a different shape from 
that of a spheroid. This could in principle introduce non axisymmetric de- 
formation in the structure of the body which might generate gravitational 
radiation. However a full significance of the result obtained would become 
clear only after a more detailed analysis which takes into account inhomo- 
geneity in the fluid configuration as well as from more realistic equations of 
state, than that of constant density used in the above analysis. 
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Figure Captions 



Fig. 1 Plots for centrifugal force F c f in units (G 2 J 2 / 'c 4 R 5 S ) for decreasing values 
of radius R in terms of Schwarzschild radius R s . 

Fig. 2 This shows two curves of ellipticity. The solid line corresponds to our 
calculation e and the dotted line represents the €h-t as used by Chan- 
drasekhar and Miller. Both the quantities are in units of (G 2 J 2 /c 6 Rj). 



Table Captions 

Table 1 Shows the ellipticity e (equation 28), e H -T (equation 29) and the cen- 
trifugal force F c f (equation 25) (units of these quantities are described 
in equation (31)) for a sequence of decreasing radius with conserved 
mass and angular momentum for a homogeneous distribution. 
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